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^ Abstract: 

^ 3+1-dimensional free inviscid fluid dynamics is shown to satisfy the criteria for exact integra- 
bility, i.e. having an infinite set of independent, conserved quantities in involution, with the 
^ ■ Hamiltonian being one of them. With (density dependent) interaction present, distinct infinite 
CN • serieses of conserved quantities in involution are discovered. Clebsch parametrization of the 
^ ' velocity field is used in the the latter analysis. Relativistic generalization of the free system is 
\Q \ also shown to be integrable. 
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Oh. The necessary conditions for a (classical) field theory to be completely integrable in more 
r-| than one spatial dimensions are expected to correspond to their 1+1-dimensional counterpart: 
j>! \ existence of an infinite number of independent conserved quantities in involution, (meaning that 
they commute at the Poisson Bracket (PB) level) with the Hamiltonian being one of the above 
These features generally induce soliton like solutions, which are physically interesting, 
both in (non-relativistic) condensed matter theory and in (relativistic) high energy physics, for 
simulating particle- like excitations. However, examples of integrable models in higher dimension 
are quite rare 0. In the present Letter, we study various aspects of integrability in the inviscid 
fluid model in 3+1-dimensions, both with and without interaction [Q]. Specifically, we show 
that (i) the free fluid models, compressible and incompressible, are integrable, (ii) compressible 
fluid model with (density dependent) interaction has two distinct infinite serieses of independent 
conserved quantities in involution and (iii) a relativistic generalization of the free fluid model is 
also integrable. "Free" model means that there is no pressure term in the equations of motion. 

Exact integrability in 1+1-dimensions has been proved [^] by exploiting the connection 
between fluid models with (Chaplygin type of) density dependent interactions and brane system 
P]. Conservation laws for irrotational fluids in higher dimensions have also been studied 

The Hamiltonian system with the basic PBs and equations of motion are the following [Q, 

1 d dv/ — dv- 

H = -pViVi + R{p) , {v,{x),p{y)} = —6{x-y) , {vi{x),Vj{y)} = — ' ^ ^ ^ ' 6{x-y). (1) 

dR 

dtp = {p{x),H} = di{pvi) , dtVi = {vi{x), H} = {vjdj)vi + (2) 
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where p and Vi denote the density and velocity fields respectively, i = 1,2,3 and a density 
dependent pressure term has been considered. Immedietely we notice that for the free theory, 
(i? = 0), for arbitrary smooth functions F{vj) and G{vj), 

dt{pFiv,)) + diipF{v,y)=0, (3) 

{ J d'x{pF),J d'yipG)} = J d'xd.ipF^) = 0. (4) 

d^) constitute an infinite number of conservation laws and (^) shows that the corresponding 
charges are in involution. Obviously the free Ti belongs to this set. Hence the free theory 
satisfies the criteria of exact integrability. However, in presence of the interaction, i? 7^ 0, the 
conservation law does not hold 0. To study the effects of interaction a better framework is 
provided by the Clebsch parametrization of the velocity field |0, Vi{x) = di9{x) + a{x)di(3{x), 
with the non-zero PBs 0, 

a 1 
{e{x),a{y)} = 6{x - y) , {f3{x),a{y)} = -6{x - y) , {e{x),p{y)} = 6{x - y). (5) 

P P 
Time evolution of the Clebsch variables are 

V- OR 

a = Vidia ; $ = Vidi(3 ; 9 = Vi{di9 - tt) + ^ ! P = 9i{pvi). (6) 

2 op 

Using above relations, one can easily check that 

dt{A{a)p) = d^{A{a)pv^) ; dt{B{f3)p) = d^{B{P)pv,) ; dt{C{a,P)p) = d^{C{a, (3)pv,), (7) 

where A, B and C are smooth functions. The above constitute the three independent serieses 
of local current conservation equations. Clearly the conserved charges A[a\ = J[A[a)p) etc. 
are trivially involutive, 

{A,[a],A2[a]} = {B,[P], B,[P]} = 0, (8) 

but the rest of the PBs between charges are non-vanishing. Hence one can take either the set 
A[a] or the set B[l3] as the independent and conserved charges in involution. This constitues 
analysis of the interacting case. Notice that Ti. does not belong to any of the three sets as it 
contains 9. It would be interesting if a bi-Hamiltonian structure could be found, in which a 
term from one of the above serieses will act as Hamiltonian to reproduce via a PB algebra 
other than (H). 

In this connection, the incompressible fiuid turns out to be integrable. We introduce the 
constraint p{x) — po ~ which generates another constraint diVi ~ due to the equations 
of motion (||). One can replace 9 = —{didi)^^dj{adjj3) from Ti. so that it belongs to C{a,(3). 
Explicitly performing the Dirac analysis |10| of the above second class pair of constraints [11] 



one can check that the conservation laws in and the corresponding charges become p- 
independent. Hence the Hamiltonian and either the set A{a) or the set B{(]) constitute the 
conserved quantities in involution. 

Some peculiarities are involved in the Noether description of the above conserved currents 
(1^) from the Lagrangian 

C = 9p + 0ap- n, (9) 
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which correctly reproduces the PBs (^) by Dirac constraint analysis or symplectic procedure 



r2| . Possibly due to the non-canonical and first order nature of the system, the action respects 
the symmetry transformations leading to the currents in a restricted way, where the arbitrary 
functions A(a), B{(3) and C(a,/3) have to be replaced by polynomials a" etc 0]. 
In the relativistic generalization, the free Lagrangian is expressed as IH 



Crei = fa, - {fj,)-^ ; a, = 0,6 + ad,(3 , f = {p, pv'). (10) 

Expanding the square root as p{l + v^Vi)^ ~ p(l + ^v^Vi + ...), and dropping the uninteresting 
/ p term [|T3|, (since it can only influence the time evolution of 6' by a constant translation), we 
can recover the non-relativistic Lagrangian. We emphasise that the symplectic structure and 
the associated PB algebra remains unaltered from (R). The Hamiltonian now is modifled to 



1 , 



Krel= p[ViVi + {l-ViVi)2], (11) 

which changes the equations of motion to the following: 

a = Lidia ; /5 = Lidi(3 ; 9 = -Liadi(3 + [viVi + (1 - ViVi)^ ; p = di{pLi), 

= = v,[2 - {1 - v^v^K (12) 

P OVi 

Notice that in the lowest order, Li ^ Vi + 0{v^), the previous equations are recovered. Clearly 
the previous conservation laws will remain intact by replacing Vi by Lj in the expressions (0) for 
the currents. Integrability is proved since in analogy to (0), we have dt{pF{vj)) + di{pF{vj)U) = 
0, (§) remains unchanged and Hrei belongs to F{vi). 

Regarding the possibility of soliton solutions in the fluid model, a comment is in order. The 
Virial theorem [|I4| proving absence of solitons in models in higher dimensions, involving only 
scalar flelds, depends on the Lagrangian having the canonical relativistic form C = d,(pid^(j)i — 
U{(j)j). The relativistic Lagrangian (|10D and its non-relativistic limit (P) both are flrst order in 
nature with a non-canonical PB algebra between the basic scalar flelds p, 9, a and p. Hence, 
a priori it is not obvious whether the Virial theorem is applicable here. A natural starting 



point is provided in the Batalin-Tyutin [T^ extension of this model (see Ghosh in §]), where 
the PBs are canonical. All the above results are derivable here in a modifled form. Also this 
version is suitable for quant izt ion. 

To conclude, we have shown that the 3+1-dimensional free fluid model satisfles the exact 
integrability criteria. For a density dependent interaction, there exist three distinct inflnite sets 
of independent, conserved charges of which terms in two sets are in involution. Subtleties in 
the Noether description of the conservation laws are pointed out. We have also shown that 
these conservation laws survive, albeit with modiflcations, in a relativistic generalization of 
the free theory. Establisment of a bi-Hamiltonian structure in case of interactions and soliton 
solutions in the model will be worth studying. How far these conclusions can be generalized to 
d-brane systems remains an open question. From another point of view, stability analysis of 



non-linear systems [0 in Clebsch representation will be interesting since, (as stated in [[16|), 
a large amount of arbitrariness in the structure of the conserved quantities, which is present 
here, is required in the analysis. 

It is a pleasure to thank Dr. B. Basu Mallick for discussions. Also I thank Dr. D. Bazeia 
for correspondence. 
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